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ABSTRACT 


This  thesis  presents  an  interpretation  of  the  Theory  of 
Biorhythms  and  develops  a  combined  statistical  model  for  a 
statistical  validation  of  the  theory.   The  model  is  then  used 
to  investigate  a  set  of  intellectual  data  and  a  set  of  acci- 
dent data.   The  intellectual  data  is  based  upon  112  academic 
grades  taken  from  four  postgraduate  school  students  over  a 
14  month  period.   The  accident  data,  taken  from  insurance 
abstracts,  consists  of  66  accidents  which  occurred  at  a  pulp 
plant  in  British  Columbia,  Canada. 
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I.   INTRODUCTION 

It  is  currently  suggested  in  scientific  theory  that  life 
functions  of  all  living  organisms  are  characterized  by  some 
form  of  cyclic  nature.   Many  of  the  cyclic  rhythm  theories 
based  on  observed  human  performance  have  been  accepted  as 
fact  while  others,  lacking  strong  scientific  and  mathematical 
proof,  have  been  treated  with  some  skepticism. 

The  most  well-known  and  best  documented  of  observed 
human  isochronic  phenomena  are  the  circadian  or  daily  rhythms. 
Circadian  rhythms  vary  in  length  from  twenty  to  twenty-eight 
hours  depending  upon  the  individual  and  environmental  cir- 
cumstances.  Included  in  these  daily  or  circadian  rhythms 
are  life  function  or  regulatory  cycles,  biochemical  cycles, 
and  socio-cultural ,  work-eat-sleep  cycles.   Extensive  re- 
search in  the  area  of  circadian  rhythms,  as  reported  by 
Squibb  Pharmaceutical  Co.  [1972],  concluded  that  regulatory 
and  biochemical  cycles  are  extremely  rigid.   Lengthy  experi- 
mentation with  circadian  work-eat-sleep  cycles  has  been  con- 
ducted by  a  variety  of  individuals  and  institutions  including 
the  Max  Planck  Institute  in  Germany  [Palmer] .   Their  research 
originated  from  two  contesting  theories.   One  theory  states 
that  the  circadian  patterned  work-eat-sleep  cycle  is  a  prod- 
uct of  socio-cultural  customs  while  the  other  theory  suggests 
it  resulted  from  the  interaction  of  the  multitude  of  circadian 
biochemical  and  regulatory  cycles  [Luce,  1973].   Scientists 
have  been  conducting  experiments  in  this  area  for  the  past 
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thirty- four  years.   The  experimental  prerequisite  to  remove 
artifical  and  natural  indicators  of  time,  such  as  clocks 
and  the  rising  and  setting  of  the  sun,  stimulated  Dr. 
Nathaniel  Kleitman  to  conduct  his  original  cave  studies,  in 
Mammouth  Cave,  Kentucky  [Luce,  1972c]  .   In  1962,  further 
studies  were  conducted  by  Michel  Siffre  who  spent  sixty- 
three  days  on  a  subterranean  glacier  in  Scarasson  Cavern  in 
the  French  Alps.   Finally,  in  1966,  Dr.  John  Mills  analyzed 
David  Lafferty's  one-hundred  and  twenty-seven  day  isolation 
in  a  cave  in  Cheddar,  England.   The  most  controlled  experi- 
ments were  conducted  at  the  Max  Planck  Institute  where 
eighty-five  subjects  lived  for  three  weeks  in  luxurious 
underground  efficiency  apartments. 

All  research  yielded  the  same  results  that  although  man 
can  artifically  control  the  length  of  his  work-eat-sleep 
cycle,  the  regulatory  and  biochemical  cycles  maintain  their 
ritualistic  twenty-four  hour  schedules.   Research  has  also 
indicated  that  the  regulatory  and  biochemical  circadian  rhythms 
influence  physiological  responses  to  environmental  stimuli. 
For  example,  the  body  reacts  differently  to  such  foreign 
agents  as  noise,  antibiotics,  hypnotics,  anaesthetics,  and 
liquor  at  different  times  of  the  day  [Squibb,  1972].   It  has 
often  been  suggested  that  circadian  rhythms  are  a  partial 
explanation  for  the  commonly  observed  circadian  post-surgical 
death  cycles  and  spontaneous  birth  cycles.   Since  circadian 
rhythms  are  constantly  observed,  the  fact  that  they  do  exist 
is  generally  accepted,  but,  in  some  cases,  why  they  exist  can 
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only  be  hypothesized  since  the  precise  roles  of  many  body 
functions  remain  a  mystery  to  the  medical  community. 

Cycle  analysis  has  not  been  limited  to  the  study  of 
circadian  rhythms.   Recently  behavioral  scientists  have  be- 
come involved  in  the  analysis  and  interpretation  of  patterns 
of  overt  human  performance.   A  considerable  amount  of  work 
has  been  conducted  by  physiologists  and  psychologists  in 
the  area  of  cyclic  phenomena  of  human  behavior.   Development 
and  application  of  several  emotional  measurement  devices  has 
led  to  various  mood  cycle  discoveries.   Tedious  experimen- 
tation and  research  by  Dr.  Max  Liischer  of  Germany  led  to  the 
original  development  of  the  Liischer  Color  Test  in  1947  which 
has  since  been  extensively  refined.   Its  validity  is  based 
upon  several  different  statistical  analyses  [Scott,  1971]. 
Adaptability,  ease  of  application,  and  lack  of  rigidity  has 
made  it  an  ideal  and  instrumental  tool  of  many  behavioral 
scientists,  from  psychologists  to  criminologists.   Although 
investigation  and  clinical  observations  by  Drs.  Halberg  and 
Panofsky  revealed  a  wide  spectrum  of  mood  cycles,  they  re- 
ported the  most  statistically  significant  periodicity  was 
seven  days  [Luce,  1972b,  Brinker,  1973].   Further  mood  cycle 
study  was  conducted  by  J.  A.  Dorland  [Brinker,  1973].   Using 
the  Dorland  Mood  Rating  Scale,  a  categorical,  subjective, 
emotional,  measurement  device,  Dorland  also  noted  the  pres- 
ence of  a  strong  seven  day  mood  cycle.   Additional  investiga- 
tion of  emotional  cycles  has  led  to  various  other  mood  scales 
such  as  Mersey's  Rating  Scale  for  Emotions  [Brinker,  1972]. 
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Additional  research  in  the  areas  of  emotional  distur- 
bances led  to  more  bizarre  theories.   Statistical  studies, 
conducted  by  Dr.  Arnold  Lieber  of  the  University  of  Miami's 
Psychiatry  Department,  based  on  1800  murders  in  Dade  County, 
Florida,  show  a  cyclic  murder  rate  coincidental  with  the 
rise  of  the  full  moon  [Hackler,  1969]  .   Doctors  from  the 
Douglas  Hospital  in  Montreal,  Canada  have  discovered  cor- 
relation between  cycles  of  intense  hostility  and  violence 
among  mental  patients  and  geomagnetic  disturbances  caused 
by  sunspots.   Many  physicians  such  as  Dr.  M.  Arborlieus  of 
Halmstead,  Sweden,  believed  that  the  awareness ,  proper  in- 
terpretation, and  application  of  physiological  and  psycho- 
logical behavior  patterns  can  be  of  significant  importance 
in  accurately  diagnosing  certain  diseases  such  as  ulcers, 
cancer,  and  psychosomatic  illness  [Luce,  1972a]. 

One  of  the  older  cyclical  theories,  based  on  observation 
of  human  behavior,  is  the  theory  of  biorhythms  developed  in 
the  late  19th  century  by  Dr.  Herman  Swoboda  of  Switzerland 
and  Dr.  Wilhelm  Fliess  of  Germany  [Thommen,  1973].   Based 
on  research  in  psychology  and  periodicity,  the  theory  states 
that  man  is  characterized  by  three  inherent  cycles  called 
biological  rhythms  or  biorhythms.   The  three  cycles,  which 
include  a  23-day  physical  biorhythm,  a  28-day  emotional 
biorhythm,  and  a  33-day  intellectual  biorhythm,  are  indica- 
tive of  the  potential  performance  capacity  of  an  individual. 
It  is  hypothesized  that  all  three  biorhythms  begin  at  birth 
and  terminate  with  death.   Swoboda  and  Fliess'  cyclic 
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hypotheses  were  based  on  pure  observation.   Since  the 
original  work  of  Swoboda  and  Fliess,  various  people  have 
explored  biorhythms  resulting  in  a  variety  of  conclusions 
and  interpretations. 

Wallenstein  and  Roberts  [1973]  tested  biorhythm  validity 
against  the  outcome  of  the  1972  UCLA-USC  football  game  along 
with  selected  historical  events  including  Custer's  prognosis 
at  Little  Big  Horn  and  Babe  Ruth's  sixty  homerun  season. 
Thommen  [1973]  attacked  validity  in  a  similar  fashion  by 
analyzing  specific  historical  events  biorhythmically  fol- 
lowed by  a  subjective  interpretation  of  the  results.   Willis 
[1972]  analyzed  data  which  included  hospital  records,  single 
car  accident  fatalities  and  athletic  events.   He  also  dis- 
cussed several  cases  of  current  application  including  the 
supposed  success  of  the  Ohmi  Railway  Company  of  Japan  in 
achieving  2,000,000  accident  free  kilometers  since  the  im- 
plementation of  biorhythms  in  driver  management.   Although 
the  chosen  specific  examples  have  corresponded  favorably  to 
biorhythmic  prediction,  this  type  of  analysis  does  not  con- 
stitute a  proper  mathematical  validation  of  the  theory. 
Closer  scrutiny  of  Ohmi ' s  accomplishments  may  not  be  so  sur- 
prising if  an  analyst  properly  applied  stochastic  processes. 
A  time  or  spatial  poisson  counting  process,  whicli  counts 
events  through  time  in  a  poisson  fashion,  could  have  been 
applied  to  predict  the  probability  of  achieving  2,000,000 
accident  free  kilometers.   Using  an  estimated  accident  rate 
based  on  past  performances  it  may  be  entirely  possible  that 
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this  could  have  occurred  with  relatively  high  probability, 
prior  to  the  use  of  biorhythms.   The  following  brief  mathe- 
matical formulation  illustrates  the  application  of  a  sto- 
chastic process  to  the  investigation  of  the  Ohmi  situation. 
To  be  complete  it  requires  the  appropriate  data,  from  the 
Nagahama  Service  Center,  needed  to  estimate  the  accident 
rate.   However, 

P(n  accidents  occur  in  the  interval  (s,s+t)  = 

-At  sy    -^n 
P(N(s  +  t)  -  N(s)  =  n)  =  P(N(t)  =  n)  =  - — n!-AtJ   . 

For  the  Ohmi  situation:  n  =  0 

t  =  1.083  years 

X   =    #  accidents/given  time  period 

P(N(1.083)  =  0)  = fXtJ   =  e-1-083A. 

Details  of  calculations,  theoretical  assumptions,  and 
the  associated  postulates  are  contained  in  "Introduction  to 
Probability  Models"  by  Ross  [1973]. 

The  absence  of  a  methodological  approach,  justification 
for  chosing  data  categories  and  acceptable  mathematical  tests 
resulting  in  statistically  significant  quantifiable  results, 
has  led  to  incomplete  results  and  subjectively  interpreted 
conclusions . 
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II.   STATEMENT  OF  PROBLEM 

The  primary  objective  of  the   thesis  was  to  establish 
the  basis  for  a  statistical  validation  of  the  Theory  of 
Biorhythms .   It  can  be  assumed  that  the  general  lack  of 
acceptance  and  understanding  of  the  theory  is  the  cause  of 
the  difficulty  to  locate  literature  on  the  subject.   The 
literature  available  is  scarce  and  the  subject  matter  is 
treated  subjectively;  consequently,  a  statistical  analysis, 
in  depth,  would  require  careful  interpretation  of  the  theory. 
The  theory  has  not  been  explicitly  stated  in  detail  in  any 
of  the  literature  while  conflicting  assumptions  supporting 
each  version  of  the  theory  depend  upon  the  author.   One 
problem  area  has  been  determining  what  a  biorhythm  actually 
measures.   Willis'  [1972]  interpretation  suggests  that  bio- 
rhythms are  indicators  of  an  individual's  relative  degree  of 
potential  performance  capacity  for  a  given  time  frame.   Other 
potential  problem  areas  include  the  interpretation  of  dif- 
ferent biorhythm  interactions  when  analyzing  observed  phe- 
nomena and  proper  classification  of  events  as  either  physical, 
emotional,  or  intellectual.   For  example,  Thommen  [1973] 
classifies  agility  as  an  intellectual  event  while  Willis  im- 
plies that  it  is  a  physical  event.   Improper  classification 
results  in  erroneous  conclusions.   It  is  often  the  case  that 
emotional  interactions  contribute  heavily  to  outcomes  of 
athletic  contests,  consequently,  athletic  contests  can  not 
be  analyzed  strictly  as  a  physical  event.   Most  of  the 
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literature  available  consists  of  short  interpretations  of 
biorhythms  and  attempts  to  draw  conclusions  about  the  theory 
based  upon  subjective  correlations  between  biorhythmic  pre- 
dictions and  special  cases  of  actual  observed  phenomena. 
The  lack  of  a  strong  statistical  foundation  is  one  major 
stumbling  block  to  universal  acceptance  of  the  Theory  of 
Biorhythms.   Statistical  validation  of  the  theory  could  pos- 
sibly lead  to  revolutionary  achievements  in  education, 
athletic  strategies,  and  industrial  safety. 
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III.   INTERPRETATION  OF  THE  THEORY  OF  BIORHYTHMS 

Human  behavior  is  characterized  by  three  inherent  cycles 
called  biological  rhythms  (biorhythms)  which  include  a  23- 
day  physical  biorhythm,  28-day  emotional  biorhythm,  and  a 
33-day  intellectual  biorhythm.   Each  biorhythm  represents 
the  individual's  relative  degree  of  potential  performance 
capacity.   It  neither  predicts  what  will  happen  in  the 
future  nor  what  happens  in  the  present.   Rather,  biorhythms 
are  indicators  of  the  individual's  physical,  emotional,  and 
intellectual  performance  capacity.   In  theory,  it  measures 
whether  the  potential  capacity  is  greater  than  or  less  than 
some  neutral  standard  on  particular  days  within  the  cycle. 
That  is,  biorhythms  measure  the  day  to  day  variation  in 
performance  potential.   The  neutral  standard  represents  the 
individual's  neutral  or  average  performance  potential  and 
is  different  for  different  individuals.   A  realization  of 
performance  potential  occurs  when  a  subject  encounters  a 
challenging  situation.   A  challenge  allows  him  to  identify 
his  full  potential.   Average  days  are  days  on  which  the  in- 
dividual does  not  encounter  challenging  situations  which 
require  a  display  of  emotional,  intellectual,  and/or  physical 
potential.   For  this  reason,  the  individual  should  be  aware 
of  the  different  catalysts  that  may  trigger  potential  re- 
actions capable  of  being  measures  by  a  biorhythm. 

Figure  1  is  a  graphical  interpretation  of  The  Theory  of 
Biorhythms . 
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Physical  Biorhythm 


Emotional  Biorhythm 


Intellectual  Biorhythm 


c 

=  Crit 
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hypo 
dina 

m  p 

I 

=  Good 

pote 
with 

I 

=  Poor 

pote 
with 

Av 

erage 

not 

can 

ical  Days  =  A  day  on  which  the  individual's  potential 
city  is  considered  unstable.   On  critical  days,  it  is 
thesized  that  individuals  can  either  have  an  extraor- 
ry  increase  in  performance  or  an  extraordinary  decrease 
erf ormance . 

Performance  Days  =  Days  on  which  there  is  a  greater 
ntial  for  above  average  performance.   These  days  occur 
in  the  first  ^  cycle. 

Performance  Days  -  Days  on  which  there  is  a  greater 
ntial  for  below  average  performance.   These  days  occur 
in  the  second  h   cycle. 

Performance  Day  =  A  day  on  which  the  individual  has 
encountered  a  challenging  situation.   Average  days 
occur  anywhere  within  the  cycle. 


Figure  1.   Graphical  Illustration  of  Biorhythms 
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Thommen  [1973]  suggests  the  following  subcategories  are 
contained  in  each  cycle. 

Physical  Cycle         Emotional  Cycle  Intellectual  Cycle 

physical  strength      sensibility     °  intelligence 

endurance  nerves  memory 

energy  feelings  mental  alertness 

resistance  intuition  logic 

confidence  cheerfulness  reasoning  power 

moodiness  reaction 

creative  ability  agility 

ambition 

Table  I.   Subcategories  of  Each  Biorhythmic  Cycle. 
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IV.   SELECTION  OF  THE  STATISTICAL  MODEL 

Although  a  variety  of  the  statistical  techniques  avail- 
able can  be  used  to  investigate  biorhythm  theory,  model  se- 
lection is  dictated  by  characteristics  of  the  data,  as  well 
as  the  appropriateness  and  quality  of  the  model. 

A.   CHARACTERISTICS  OF  THE  DATA 

Three  approaches  are  available.   The  first  approach 
consists  of  collecting  data  resulting  from  a  phenomena  which 
is  predictable  from  biorhythm  theory  and  then  selecting  a 
statistical  technique  whose  assumptions  are  satisfied  by 
the  data  and  whose  results  most  appropriately  illustrate  the 
validity  of  the  theory  of  biorhythms.   In  using  this  approach 
it  should  be  obvious  that  the  choice  of  the  technique  is  in- 
fluenced by  the  characteristics  of  the  data,  which  must  con- 
form to  the  assumptions  of  the  statistical  model.   For 
example,  the  most  common  assumption  of  parametric  tests  is 
homogeneity  of  sampling  units.   Although  it  may  be  possible 
to  restore  homogeneity  through  transgeneration ,  the  use  of 
non-homogenous  data  in  a  parametric  test  leads  to  invalid 
results.   Consequently,  the  absence  of  proper  data,  which 
can  force  the  use  of  a  less  appropriate  statistical  test, 
may  result  in  inconclusive  evidence  supporting  the  theory. 

The  second  approach  consists  of  selecting  the  most  ap- 
propriate statistical  tests  which  best  illustrate  the  validity 
of  the  biorhythm  theory,  then  collect  data  conforming  to  the 

assumptions  of  that  test. 
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A  third  approach  was  considered  out  of  necessity.   Data 
was  often  difficult  to  obtain  which  both  illustrated  par- 
ticular aspects  of  the  theory  and  conformed  to  the  assump- 
tions of  the  chosen  statistical  model.   Consequently,  it 
was  necessary  to  transgenerate  raw  data  into  categorical 
variables  and  then  proceed  with  a  non-parametric  test. 

B.   APPROPRIATENESS  OF  STATISTICAL  MODEL 

Appropriateness  of  the  statistical  test  is  determined 
by  the  objective  of  the  thesis.   Various  parametric  and 
non-parametric  tests  were  considered.   The  robust  t-test 
analysis  of  the  variance,  and  their  non-parametric  counter- 
parts would  have  been  inappropriate  since  they  allowed  anal- 
ysis of  only  specific  segments  of  the  cycle.   Consequently, 
even  if  the  results  were  favorable,  the  evidence  would  be 
inconclusive.   For  example,  depending  upon  homogeneity 
qualities  of  the  data,  analysis  of  the  variance  could  have 
been  used  to  compare  the  statistical  differences  between 
biorhythmically  predictable  critical  days  and  the  actually 
observed  critical  days.   However,  this  technique  would  not 
allow  a  simultaneous  investigation  of  the  statistical  sig- 
nificance of  the  occurrence  of  other  categorical  events  be- 
tween the  biorhythmic  critical  days.   The  results  of  the 
statistical  models  must  completely  support  the  objective 
of  the  thesis.   Therefore,  tests  which  offer  only  a  partial 
analysis  were  considered  inappropriate. 

An  appropriate  statistical  model  for  testing  the  validity 
of  biorhythms  should  test  the  existence  of  a  relationship  or 
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dependency  between  observed  good  performance  days  and  bio- 
rhythmically  predictable  good  performance  days  and  observed 
poor  performance  days  and  biorhythmically  predictable  poor 
performance  days.   Figure  2  illustrates  this  relationship. 

It  is  intuitively  appealling  that  a  relationship  exists 
since  most  of  the  good  performance  days  occur  during  the 
first  %  cycle  and  most  poor  performance  days. occur  during 
the  second  %  cycle.   An  appropriate  statistical  model  should 
be  able  to  determine  the  presence  of  a  relationship  or  de- 
pendency between  theoretically  predictable  days  and  observed 
days  as  illustrated.   If  the  occurrence  of  good  days  and 
poor  days  are  equally  likely  to  occur  anywhere  within  the 
entire  cycle  then  the  distribution  of  these  days  should  be 
uniform.   Figure  3  illustrates  this.   A  second  statistical 
test  should  illustrate  uniformity  or  non-uniformity.   If  the 
distribution  is  not  uniform,  this  test  should  be  able  to 
identify  the  distributional  form.   Consequently,  another 
statistical  model  should  be  used  to  illustrate  uniformity 
or  non-uniformity. 

C.   QUALITY  OF  STATISTICAL  MODEL 

When  alternative  statistical  models  are  equally  appro- 
priate, the  technique  with  greater  quality  is  chosen.   Al- 
though there  may  be  various  ways  to  determine  the  quality 
of  a  model,  the  most  objective  measure  of  quality  is  power- 
efficiency.   Suppose  two  tests,  A  and  B  are  equally  appro- 
priate.  If  test  A  requires  N   samples  to  have  the  same  power 
as  test  B  has  with  K,     samples  then  test  B  has 
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A  =  Theoretical  envelope  of  biorhythmically  predictable  good 
days  ; 

B  =  Hypothetical  distribution  of  observed  good  days; 

C  =  Theoretical  envelope  of  biorhythmically  predictable  poor 
days  ; 

D  =  Hypothetical  distribution  of  observed  poor  days. 


Figure  2.   Theoretical  Biorhythm  Superimposed  Upon  Hypothetical 
Distribution  of  Observed  Good  and  Poor  Performance 
Days. 


POWER-EFFICIENCY  =  (100) 


ft) 


where  power  is  defined  as  the  probability  of  rejecting  the 
null  hypothesis  when  it  is  false  [Siegel,  1956]. 

Considering  the  characteristics  of  the  data  available, 
and  appropriateness  and  quality  of  the  statistical  models, 
the  Chi -Square  Contingency  Test  for  independence  of  criteria 
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A  =  Uniform  distribution  of  observed  good  days; 
B  =  Uniform  distribution  of  observed  poor  days. 

Figure  3.   Hypothetical  Sample  from  a  Uniform  Distribution 
of  Good  and  Poor  Performance  Days. 

and  the  Chi-Square  Goodness -of -  Fit-Test  for  distributional 
form  were  chosen  as  primary  statistical  techniques.   Secon- 
dary techniques  are  available  to  give  supporting  evidence. 
When  sample  size  is  so  small  that  the  Chi-Square  Contingency 
Test's  expected  value  criterion  does  not  hold,  the  Fischer 
Test  for  exact  probability  (without  Tocher's  modification) 
should  be  used.   When  the  Chi-Square  Goodness -of -Fi t -Test ' s 
expected  value  criterion  does  not  hold  and  category  combina- 
tion reduces  the  degrees  of  freedom  to  an  undesirable  level 
the  Kolmorgorov-Smirnov  Test  for  distributional  form  should 
be  used. 
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These  statistical  techniques  indicate  the  presence  or 
lack  of  presence  of  a  relationship  and  define  the  distribu- 
tional form.   They  do  not  attempt  to  explain  why  the  rela- 
tionship is  or  is  not  present.   Elaborate  linear  and 
curvilinear  regression  techniques  are  available  for  the  de- 
velopment of  quantitative  explanatory  models. 

In  addition  to  analyzing  the  data  with  r.espect  to  de- 
pendency and  distributional  form,  it  would  also  be  desirable 
to  measure  the  degree  of  correlation  between  the  theoretical 
biorhythm  and  the  actual  observations.   This  can  be  most 
easily  accomplished  by  utilizing  the  Chi-Square  Contingency 
Coefficient. 

D.   CHI -SQUARE  CONTINGENCY  TEST 
1.   Theory 

A  two-way  Chi-Square  Contingency  Test  was  used  to 
illustrate  the  dependency  between  the  day  of  occurrence  of 
the  observed  event  and  the  sequential  day  of  the  cycle. 
Biorhythm  theory  states  that  an  individual's  capacity  to 
perform  is  dependent  upon  the  sequential  day  of  the  cycle. 
Specifically,  there  is  a  greater  capacity  for  an  individual 
to  excel  during  the  first  half  cycle  and  a  reduction  in 
capacity  during  the  last  half  of  the  cycle.   Therefore,  the 
Chi-Square  Contingency  Test  was  used  to  investigate  the  de- 
pendency between  the  day  of  occurrence  of  observed  human 
performance  and  the  sequential  day  of  the  cycle. 

a.   Data  is  arranged  in  the  following  tableau: 
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criteria  1 
12  3  4  5.  . 


1 

2 

3 

criteria  4 

2     5 


ij 


re 


'I 


b.  Given  a  sample  of  n  observations. 

c.  Each  observation  can  be  classified  according  to 
two  different  criteria.   One  criteria  has  c  levels  with 

c  >  2  and  the  other  criteria  has  r  levels  with  r  >  2. 

d.  c    r 

E    E   x. .  =  n. 
j=l  i=l   ^ 


e.  c    r 

E    E   P. .  «  1. 
j  =  l  i-1   ^ 

f.  nP.  .  >  5. 

g.  Let  x.   be  the  number  of  observations  which  fall 
into  class  i  criteria  one  and  class  j  criteria  two. 

h.   Let  P.   be  the  probability  that  an  observation 
falls  into  class  i  criteria  one  and  class  j  criteria  two. 
i .   .       x. 


P. 
i 


where  x. 

n  i 


c 

E   x. 


j  =  l 


x  . 

— '-2-     where  x 
n  .  3 


r 

E   x. 
i-1   1J 
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k.   H   :  P. .  =  (P.  ) (P  .)  for  all  i  and  j 


1. 


o 
H,  :  P.,  f    (P,  )(P  ,) 


J 


If  H  is  true : 
o 


c    r 
V  =   Z    £ 
j=l  i=l 


(x. .  -  nP.  P  .) 
nP.  P  . 


=  X2(r-l)(c-l) 


(1-a) 


m.   Reject  H   if: 
J  o 


V  >  X2(r-l)(c-l) 


(1-a) 


Procedure  for  the  Application  to  the  Theory  of 
Biorhy thms 

a.   Data  tableau 


Criteria  1 


1st  H 

Cycle 

Criteria 

2      A 
2nd  h 

Cycle 


Observed 
Good  Days 


Observed 
Bad  Days 


xll 

X12 

X21 

X22 

1  h   cycle: 

_nd  ,     , 

2  h   cycle: 


represents  biorhythmically  predictable 
good  performance  days 

represents  biorhythmically  predictable 
poor  performance  days . 


b.   x.   represents  the  number  of  observations  which 
are  categorized  by  both  criteria  levels  i  and  j.   For  example, 
x?1  represents  the  number  of  observed  good  performance  days 
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on  which  the  individual's  performance  was  above  average, 
which  occurred  during  the  second  half  of  the  cycle. 

Figure  4  is  a  graphical  illustration  of  x..,  xi?> 

c.   P.   and  P    represent  the  maximum  likelihood  esti- 
l  •       •  j    r 

mator  for  the  probability  that  an  observation  is  character- 
ized by  the  respective  criteria  level.   For  example,  P.. 
represents  the  probability  that  a  good  or  bad  day  occurs 
during  the  1st  half  of  the  cycle. 

d-  2  2 

Ex..  Ex.. 

x.      .  .,   11  •       .  ,   11 

P.   =  _JLL  =  3  =  1    J  p    =  i=l    J 


i  •     n      2    2  •  j     2    2 

E    E   x.  .  E    E   x.  . 

j  =  l  i=l   1J  j  =  l  i=l   XJ 

p    =        XH  +  X12  J,  Xll  +  X2i 


1-    x1]L  +  x12  +  x21  +  x22    -1   x1;L  +  x12  +  x21  +  x22 
p    =  X21  +  X22 p  X12  +  X22 


2-    xn  ♦  x12  +  x21  +  x22    -2    xn  +  x12  +  x21  +  x22 

For  best  results 

Efx.  -)    =    nP.    P    .    >    5. 
ij  i  •    •  j    - 

If   H      is    true: 
o 

2        2    (x.  .    -    nP.    P    .)2  M      ■, 

V=      E         jilL, .x'    'J         =    x2(r-l)(c-l)(1"a). 


j=l   i=l  nP.    P    . 


g.       Reiect   H      if: 
°  -  o 


,2^_-,w_^(l-a) 


V   >    X'(r-l)(c-l) 
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Theoretical  Biorhythm  superimposed  upon  hypothetical  dis- 
tributions of  observed  good  performance  and  poor  perform- 
ance days  . 


1st  %  cycle 

2nd 
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total  number  of  observed  good  performance  days  occurring 
during  the  1st  Jj  cycle.   They  appear  in  quadrant  II. 

total  number  of  observed  poor  performance  days  occurring 
during  the  1st  h   cycle.   They  appear  in  quadrant  III. 

total  number  of  observed  good  performance  days  occurring 
during  the  2nd  %  cycle.   They  appear  in  quadrant  I. 

total  number  of  observed  poor  performance  days  occurring 
during  the  2nd  %  cycle.   They  appear  in  quadrant  IV. 


Figure  4.   Location  of  Statistical  Data  Within  a  Biorhythm. 

E.   CHI -SQUARE  COEFFICIENT  OF  CONTINGENCY 

The  contingency  coefficient  is  especially  useful  in  that 
it  measures  the  degree  or  strength  of  the  relationship  which 
was  detected  using  the  Chi-Square  Contingency  Test  for  de- 
pendence of  criteria.   It  is  calculated  as  follows. 
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X2  ♦  N 


C  =  Coefficient  of  Contingency. 
X2  =  Test  statistic  value  (V)  calculated  from 
the  Chi-Square  Contingency  Test  (previously  discussed). 

N  =  Total  number  of  observations. 


F.   CHI-SQUARE  GOODNESS-OF-FIT  TEST 

If  the  occurrence  of  a  categorical  event  was  equally 
likely  to  occur  anywhere  within  the  cycle  then  the  distribu- 
tion of  these  events  would  be  uniform.   Consequently,  valida- 
tion of  the  theory  of  biorhythms  necessitates  illustrating 
that  the  occurrence  of  categorical  events  is  not  uniformly 
distributed.   Investigating  distributional  form  is  most 
easily  accomplished  by  utilizing  the  Chi-Square  Goodness-of- 
Fit  Test. 

1.   The  data  is  organized  in  the  following  tableau. 


Interval 

Freque 

ncy 

P. 

1 

nP. 
1 

a  <  x  <  b 

Yi 

pi 

EOy 

b  <  x  <  c 

Y2 

• 

P2 

• 

E(Y2) 

g  <  x  <  h 

• 

Yk 

• 

• 

E(Yk) 

2.  n:   represents  the  sample  size. 

3.  a  <  x  <  b:   represents  an  interval  of  time. 

4.  Y. :   represents  the  frequency  of  occurrence  of  the 
observed  phenomena  within  the  i    interval. 
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5.  P.:   represents  the  probability  that  an  observed 
event  will  occur,  under  the  null  hypothesis,  during  the 
associated  interval. 

6.  nP.:   represents  the  expected  number  of  observations 
within  the  i    interval  if  the  distribution  of  the  frequency 
of  occurrence  of  the  phenomena  conformed  to  the  null 
hypothesis. 

7.  For  most  significant  results  nP.  >  5. 

8.  H   :   Hypothesized  distributional  form 

0  '  r 

H1  :   The  distribution  is  other  than  as  hypothesized. 

9.  If  H   is  true : 

o 

k   (Y.  -  nP.)2  m   -> 

V-   Z     \Pm      X        =  x2Ck-r-l)<>0. 
i=l      i 

k-r-1  =  degrees  of  freedom  for  the  test  statistic 

k  =  number  of  intervals 

r  =  number  of  estimators  used  in  the  hypothe- 
sized distribution. 

10.  Reject  H   if: 

J      o 

V  >  X2(k-r-l)(1"a). 

11.  In  an  application  to  biorhythms,  the  initial  test  is 

for  the  presence  of  a  uniform  distribution. 

H   :   the  distribution  is  uniform. 
o 

H..  :   the  distribution  is  other  than  uniform 

P   -  b  -  a 

1  "  h  -  a 

P   =  d  -  c 
*2    h  -  a 


etc 
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Acceptance  of  the  Contingency  Test's  null  hypothesis 
indicates  that  the  day  of  occurrence  of  observed  events  is 
independent  of  the  sequential  day  of  the  cycle.   In  other 
words,- there  is  no  relation  between  the  days  of  occurrence 
of  biorhythmically  predictable  events  and  the  days  of  occur- 
rence of  the  observed  events.   This  is  an  indication  that 
the  distribution  of  the  occurrence  of  sampled  events  should 
be  uniform.   It  is  considered  conclusive  evidence  (statis- 
tically) that  the  Theory  of  Biorhythms  is  invalid,  if  for 
the  same  data,  H   for  the  Chi-Square  Contingency  Test  and 
the  Chi-Square  Goodness -of -Fit  Test  is  accepted.   Rejection 
of  both  null  hypotheses  suggests  validity  of  the  theory. 
Further  investigation  of  the  distributional  form  would  be 
desired.   The  presence  of  a  known  distributional  form  such 
as  a  normal  or  poisson  would  be  a  useful  predictive  device. 

G.   ESTIMATION  OF  PARAMETERS 

Investigating  the  presence  of  more  complex  distributional 
forms  usually  results  in  estimating  parameters  of  the  hypothe- 
sized distribution.   Consequently,  estimators  with  good 
properties  are  important.   The  method  of  moments  and  the 
method  of  maximum  likelihood  techniques  were  considered. 

1 .   Methods  of  Moments 

a.  Given  a  population  random  variable  X  is  dependent 
upon  k  unknown  parameters  Q-,  ,  82,...  ,0,  . 

b.  In  general  the  first  k  population  moments  will 
depend  upon  these  parameters  (i.e.,  E(X)  =  g  (G .  ,  0~  , .  .  .  ,  QjJ)  • 
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c.  Equate  the  first  k  sample  moments  M- ,U? , . . . ,M,  , 
to  the  first  k  population  moments 

1   n      k 
k   n  i=1   iJ 

d.  Solve  for  the  desired  parameter. 
2 .   Method  of  Maximum  Likelihood 

a.  Develop  a  likelihood  function 

n 

L(e1,e2,...,ek)  =  _z  fx(x). 

n 

b.  Since  .|..f  (x)  is  just  another  way  of  writing  a 

sample's  joint  distribution  function,  the  value  of  f(x)  is 
actually  a  probability. 

c.  The  idea  of  the  MLE  technique  is  to  find  the 
value  of  the  unknown  parameters  within  the  joint  density 
which  will  maximize  it.   This  value  maximizes  the  probability 
of  occurrence.   The  procedure  is  as  follows: 

n 
L(6)  =   E   f  (x) 
i=l   x 

K(6)  =  ln(L(9)) 

d(K(6))  _  d(ln(L(9)))  _ 
dG         d~6         U 

solve  for  0  and  label  it  0. 

In  most  cases  the  method  of  moments  and  the  .^thod 
of  maximum  likelihood  yield  the  same  results.   However,  in 
those  cases  for  which  the  parameter  estimators  are  different, 
the  properties  of  the  estimators  should  determine  which 
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estimator  will  be  used.   In  many  cases  it  is  desirable  to 
have  unbiased,  consistent,  estimators  with  the  least  vari- 
ance and  mean  square  error  (MSE) . 
.  3.   Bias 

is  an  unbiased  estimator  for  0  iff  E  (©)  =  o.   The 
bias  in     is  given  by  h~   =  E(®  )  -  6. 

4 .  Consistency 

is  a  consistent  estimator  for  8 

iff  limit  EC© )  =  6 

n  -*-  °° 

limit  V(®)  =  0. 
n  ■*■  °° 

5.  Cramer  Rao  Lower  Bounds 


a2   = 


nEfCgf  log  fx(x))2] 


MSE  -    CI  t   b'(6))2 _ 

nE[(g|  log  fx(x))2] 

Table  II  summarizes  the  combined  statistical  model 
suggested  for  testing  the  presence  of  a  biorhythm. 

H.   CRITICALITY 

Critical  days  are  defined  to  be  days  on  which  an  indi- 
vidual's potential  capacity  for  performance  is  unstable.  If 
an  individual  encounters  a  challenging  situation  on  a  bic- 
rhythmically  critical  day  then  it  is  possible  that  he  will 
either  have  extreme  negative  potential  or  extreme  positive 
potential.   For  example,  it  has  been  suggested  that  heart 
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Statistical  Model 


Application 


1.   x2  Contingency  Test 

(Use  Fischer  Test 
for  small  samples) 


Tests  for  the  presence  of  a 
relationship  between  the 
different  levels  of  perform- 
ance and  the  days  on  which 
they  occurred. 


2.   x2  Coefficient  of 
Contingency 


Tests  the  strength  of  the 
relationship  detected  in 
step  #1. 


X2  Goodness-of -Fit  Test 

(Use  Kolmogornov- 
Smirnov  Test  for 
small  samples) 


1st  Application: 

Test  for  uniformity  (test 
for  the  presence  of  cyclic 
phenomena) . 

2nd  Application: 

Test  for  hypothesized  dis- 
tributional form.   This 
allows  the  analyst  to  in- 
vestigate the  following: 

a)  Form  of  cyclic  phenomena 

b)  Periodicity 

c)  Frequency 

d)  Criticalities 


Note: 

Tests  should  be  applied  in  the  above  order. 
Detailed  discussion  and  application  procedures  for  the 
Fischer  Exact  probability  Test  and  the  Kolmogornov-Smirnov 
Test  for  distributional  form  can  be  found  in  Siegel  [1956]. 


Table  II.   Combined  Statistical  Model  for  Testing  the  Presence 
of  a  Biorhythm. 
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attack  victims  encountered  their  attacks  on  critical  days. 
It  is  obvious  that  a  heart  attack  can  be  classified  as  a 
critical  event.   However,  did  a  statistically  significant 
number  of  heart  attack  victims  experience  their  heart  at- 
tacks on  biorhythmic  critical  days?   Thus,  criticality  be- 
comes an  important  factor  when  investigating  the  theory. 
For  this  reason,  probabilistic  arguments  have  been  formu- 
lated for  use  in  the  following  statistical  model  presented 
for  investigating  the  significance  of  an  event  occurring  on 
a  critical  day. 

1.   Definitions 

Single  Criticality  =  day  on  which  only  one  biorhythm 
crosses  the  time  axis. 

Double  Criticality  =  days  on  which  any  two  biorhythms 
intersect  while  crossing  the  time  axis. 

Triple  Criticality  =  day  on  which  all  three  biorhythms 
intersect  while  crossing  the  time  axis. 

P(PCD)  =  probability  that  a  physical  critical  day 
occurs  by  chance. 

P(ECD)  =  probability  that  an  emotional  critical  day 
occurs  by  chance. 

P(ICD)  =  probability  that  an  intellectual  critical 
day  occurs  by  chance. 

P(PCDxECD)  =  probability  that  a  double  critical  day 
occurs  by  chance  where  the  intersecting  biorhythms  are  physi- 
cal and  emotional. 
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P(PCDxICD)  =  probability  that  a  double  critical  day 
occurs  by  chance  where  the  intersecting  biorhythms  are  physi- 
cal and  intellectual. 

P(ECDxICD)  =  probability  that  a  double  critical  day 
occurs  by  chance  where  the  intersecting  biorhythms  are  emo- 
tional and  intellectual. 

P(PCDxECDxICD)  =  probability  that  a. triple  critical 
day  occurs  by  chance. 

It  is  possible  that  there  are  differing  degrees  of 
double  criticalities .   For  example,  it  may  be  that  a  physi- 
cal x  emotional  double  critical  day  possesses  a  higher 
degree  of  criticality  than  does  a  physical  x  intellectual 
double  critical  day. 

2 .   Occurrence  of  Critical  Days 

The  following  table  illustrates  the  frequency  of 
occurrence,  in  days  of  all  types  of  critical  days. 

Criticality   Biorhythm         Frequency  (days) 

Single  Physical  (P)  11.50 

Single  Emotional  (E)  14.00 

Single  Intellectual  (I)  16.50 

Double  P  x  E  161.00 

Double  P  x  I  189. 75 

Double  E  x  I  231.00 

Triple  P  x  E  x  I  2656. 50 

Table  III.   Frequency  of  Occurrence  of  Criticalities. 

The  frequency  (days)  is  based  upon  the  respective 
cycle  length.   For  example,  a  single  criticality  for  a  physi- 
cal biorhythm  occurs  every  11%  days.   Since  the  physical 
biorhythm  is  hypothesized  to  be  a  23-day  sinusoid,  it  crosses 
the  axis  every  11%  days  (i.e.,  23/2  =  11%). 
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If  the  events  are  equally  likely  (uniform)  to  occur 
at  any  point  in  time,  then  the  following  probabilities  mea- 
sure the  likelihood  that  the  observed  event  will  occur  on 
the  given  biorhythmic  critical  day.   Calculation  was  based 
on  the  frequency  of  occurrence  of  critical  days  and  con- 
sideration was  given  to  the  uncertainity  of  birth  data. 


0  Hr     ±  6  Hr    ±12  Hr 
Chance  Probability  Variance   Variance   Variance 


P(PCD) 

.0869 

.1304 

.1739 

P(ECD) 

.0714 

.1071 

.1428 

P(ICD) 

.0606 

.0909 

.1212 

P(PCD  x 

ECD) 

.0062 

.0093 

.0124 

P(PCD  x 

ICD) 

.0053 

.0078 

.0105 

P(ECD  x 

ICD) 

.0043 

.0065 

.0086 

P(PCD  x 

ECD  x 

ICD) 

.00037 

.0005 

.0007 

Table  IV.   Chance  Probability  of  Occurrence  of  Criticalities 

3.  Selection  of  Statistical  Model  for  Criticality 
Biorhythm  data  is  naturally  dichotomous  into  those 

days  on  which  critical  events  occurred  and  those  days  on 
which  critical  events  did  not  occur.   Since  an  important 
question  to  be  answered  in  validating  the  concept  of  criti- 
cality would  be;  "Are  the  number  of  observed  critical  days 
statistically  significantly  greater  than  the  number  that 
would  be  expected  if  critical  event  days  were  to  have  oc- 
curred by  pure  chance?";  the   binomial  model  was  chosen. 

4 .  Binomial  Model 

Assume  X~Binomial  (n,p) 
H 


H. 


a 


P  =  *P(critical  day  occurs  by  chance) 
P  >  P(critical  day  occurs  by  chance) 
Level  of  significance 
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*Select  appropriate  probability  from  chance  prob- 
ability table. 

Suppose  that  x  observations  out  of  n  possible  ob- 
servations are  classified  as  critical  events. 

Define  P  (X  <  x)  =  probability  that  at  most  x  criti- 
cal events  could  have  occurred  by  chance. 

n 
P(X  <  x)  =   E  (n)  P^l-P)11"1 
i=0  x 

Decision:   If  P  <  a  reject  H   in  favor  of  H..  . 

J  o  1 

The  decision  to  reject  H   says  that  the  number  of 
critical  events  which  have  occurred  is  significant. 

For  n  >  25  it  is  suggested  that  a  normal  approxima- 
tion to  the  above  probability  be  used. 

P(X  <  x)  =  *(%—===:)• 
/nP(l-P) 

Additional  information: 

E(X)  =  expected  number  of  critical  events  under 

H   =  nP. 
o 

V(X)  =  variance  of  critical  events  under  H   =  nP(l-P). 

The  Binomial  Model  is  primarily  used  when  investigat- 
ing the  significance  of  single  criticalities .   The  Multinomial 
Model  should  be  used  when  investigating  the  statistical  sig- 
nificance of  data  which  contains  single,  double,  and/or  triple 
criticalities . 

5 .   Multinomial  Model 

Suppose  that  N  observed  events  are  partitioned  into  S 
mutually  exclusive  events  A-., A-,..., A  with'  respective 
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probabilities  of  occurrence  P..,P2,...,P  such  that  A..  + 

A-  +  . . .  +  A   =  N  and  P-  +  P_  +  . . .  +  P   =  1.   Then  P(A, 
2  s  12  s  v  1 

occurs  k.,  times;  A~  occurs  k_  times:  ..  .;  A  occurs  k 
1      '2         2  s         s 

times)  is  equal  to 

|  P*l  P^2  P^3  ...  pks 

k   k   k      k     X    2    3       s 

K1)K2,K3'' • ' '  s/ 

I.   BIORHYTHM  ALGORITHM 

1.  Note  the  date  (month-day-year)  of  the  observed  event. 

2.  Calculate  the  subject's  age  in  days,  from  the  sub- 
ject's birth  date  up  to  and  including  the  date  of  the  ob- 
served event.   The  calculation  must  include  regular  leap 
year  and  centurial  leap  year  considerations. 

3.  Let 

A  =  subject's  age  in  days, 

CL  =  cycle  length  =    23  days  for  physical  cycle 

28  days  for  emotional  cycle 
33  days  for  intellectual  cycle 

CC  =  whole  number  of  completed  cycles 

(whole  number  associated  with  A/CL) 

SD  =  the  day  of  the  current  biorhythm  cycle  on  which 
the  observed  event  occurred. 

SD  =  C  cT  "  CC) (CL) • 

J.   DATA  COLLECTION 

Statistical  validation  of  the  theory  necessitates  per- 
forming two  complete  statistical  analyses.   The  first  analysis 
should  use  empirical  data  collected  in  a  laboratory  environ- 
ment, for  each  subcategory.   Experimental  design  should  screen 
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out  as  much  exterior  influences  and  as  much  interaction  as 
possible,  from  the  subcategorical  data  being  collected. 
The  second  analysis  should  utilize  actual  field  data.   It 
is  here  that  cycle  interaction  will  have  to  be  assessed. 
For  example,  if  the  analyst  chooses  to  gather  bowling  data 
for  physical  field  data,  he  should  be  aware  that  there 
exists  some  degree  of  emotional  and  intellectual  interac- 
tion.  If  the  analyst  chooses  to  analyze  academic  test 
scores,  he  should  be  aware  of  which  tests  challenged  the 
subject's  intellect,  and  that  the  importance  of  the  test  and 
the  length  of  the  test  (i.e.,  time)  will  cause  emotional 
and  physical  interaction  respectively. 

Since  human  behavior  is  extremely  complex,  there  are  a 
multitude  of  influencing  factors.   Therefore,  it  should  be 
understood  that  biorhythms  are  not  'cure-alls'.   A  biorhythm 
is  hypothesized  to  be  merely  one  very  small  indicator  of 
the  variation  of  performance  potential.   The  many  influenc- 
ing factors  can  easily  contaminate  biorhythmic  analysis. 
Consequently,  it  is  not  expected  that  field  data  will  yield 
as  strong  a  statistical  significance  as  experimental  data. 
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V.   ANALYSIS  OF  DATA 

A.   ANALYSIS  OF  INTELLECTUAL  DATA 

Four  graduate  students  at  the  Naval  Postgraduate  School 
were  chosen  as  subjects  for  investigation  of  intellectual 
biorhythms.   Academic  grades  from  15  courses  measured  over 
a  14  month  period  provided  115  data  points  for  analysis. 
Each  student,  pursuing  his  Master's  Degree  in  Operations 
Research,  was  evaluated  with  respect  to  the  same  courses. 
To  eliminate  a  degree  of  subjectivity,  only  mathematically 
oriented  courses  were  used.   It  was  believed  that  grades 
were  arrived  at  more  objectively  and  that  the  subject's 
identification  of  his  own  neutral,  average  performance 
would  be  more  accurate.   Identification  of  this  average  per- 
formance is  considered  critical  since  it  allows  the  analyst 
to  categorize  those  performances  which  were  above  and  those 
performances  which  were  below  average.   Table  V  summarizes 
the  data  base. 

The  data  was  next  arranged  in  the  Intellectual  Data 
Tableau  which  is  shown  as  Table  VI. 

1 .   Tests  for  the  Presence  of  an  Intellectual  Biorhythm 

a.   x2  Contingency  Test 

II  :   Observed  good  performance  days  and  observed 

poor  performance  days  are  independent  of  where  they  occur 

within  the  intellectual  biorhythm. 

rL  :   Observed  good  performance  days  and  observed 

poor  performance  days  are  dependent  of  where  they  occur  within 

the  intellectual  biorhythm. 
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Avg 

Slightly 

Well 

Slig 

htly 

Well 

Cycle 

Per 

Above 

Above 

Bel 

ow 

Below 

Day 

0 

Avg  (+) 

Avg  (++) 

Avg 

CO 

Avg  (--) 

Critical  0 

1 

1 

1 

2 

2 

1 

3 

3 

1 

1 

4 

2 

5 

3 

1 

1 

6 

2 

1 

7 

1 

1 

1 

1 

8 

3 

3 

9 

2 

1 

10 

1 

1 

11 

4 

2 

12 

4 

2 

13 

9 

u 

1 

1 

14 

1 

1 

15 

1 

1 

1 

16 

2 

1 

17 

18 

1 

1 

1 

1 

19 

2 

20 

1 

3 

3 

21 

1 

3 

22 

1 

23 

2 

1 

., 

2 

24 

1 

2 

25 

2 

1 

26 

2 

2 

27 

1 

3 

1 

28 

1 

1 

1 

29 

2 

1 

30 

1 

1 

2 

1 

31 

2 

1 

1 

32 

1 

2 

33 

The  data  points  within  the  tableau  represent  the  frequency 
of  occurrence.  For  example,  there  were  four  above  average 
performances  on  the  11th  day  of  the  cycle. 


The  data  can  now  be  easily  categorized  for  use 
various  statistical  models. 


in  the 


Table  VI.   Intellectual  Data  Tableau 
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a  :   Level  of  Significance  =  .005 


Observed  Good 
Performance 
Days 


Observed  Poor 
Performance 
Days 


1st  h  Cycle 
2nd  h  Cycle 
3rd  h  Cycle 
4th  \   Cycle 


21 

4 

15 

9 

8 

14 

2 

16 

X3C.995)  =  12.8    V  =  25.7 


Table  VII.   Intellectual  Data  x2  Contingency  Test. 

Decision:   Reject  H   in  favor  of  H- .   This  in- 
J  o  1 

dicates  the  presence  of  a  relationship  between  the  different 
levels  of  performance  and  the  biorhythmic  days  on  which  they 
occurred. 

b.   x2  Coefficient  of  Contingency 


C  = 


^V  +  N     \ 


25.7 
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\  2  5.7  +  89 

This  represents  the  strength  of  the  relation- 
ship.  It  is  important  to  realize  that  C  =  .48  indicates 
that  the  relationship  is  strong.   Notice  that  even  if  the 
relationship  were  perfect,  it  would  be  impossible  for  C  =  1, 
since  V/V  +  N  will  always  be  less  than  1. 

2 .   Distribution  of  Average  Performance  Days 
a.   Test  For  Uniformity 

X2  Goodness  -  of -Pit  Test 

H  :   The  distribution  of  average  performance  days 
is  uniformly  distributed  over  the  intellectual  biorhythm. 
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H1 :   The  distribution  of  average  performance 
days  is  not  uniformly  distributed  over  the  intellectual 
biorhythm. 

a  :   Level  of  Significance  =  .01. 


Interval  (days) 

Fre 

que 

acy 

P. 

l 

nP. 

l 

1  <  t  <   8 

8  <  t  <  16 

16  <  t  <  26 

26  <  t  <  33 

5 
4 
8 
8 

8/33 

8/33 

10/33 

7/33 

6.06 
6.06 

7.57 

5.30 

1 

X|(.99)  =  11.3     V  =  2.2! 


Table  VIII.   Intellectual  Data  (Average  Performance  Days): 
X2  Goodness-of -Fit  Test  for  Uniformity. 

Decision:   Do  not  reject  H  . 

J  o 

Kolmogorov-Smirnov  Test 

Utilizing  the  same  H  ,  FL  ,  and  a  as  above: 
6  o '   1 ' 


Interval 

1  <  t  < 

Frequency 

5 

F  (x) 

8/33 

S25(x) 

5/25 

lFo  "  S25l 

.0424 

8  <  t  <  16    16  <  t  <  26   26  <  t  <  33 
4  8  8 

16/33         26/33        33/33 


9/25 


1248 


17/25 


1079 


25/25 


00.00 


Critical  D  =  .32   Calculated  D  ■  .1248. 

Tabic  IX.   Intellectual  Data  (Average  Performance  Days): 
Kolmogorov-Smirnov  Test  for  Uniformity. 

Decision:   Do  not  reject  H  . 

J  o 

Since  the  distribution  is  uniformly  distributed, 
there  is  no  need  to  further  investigate  distributional  form. 
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The  above  analysis  confirms  the  assumption  that  average  per- 
formance days  can  occur  anywhere  within  the  cycle  with  equal 
probability  (uniformity) . 

3.   Distribution  of  Above  Average  Performance  Days 

For  stronger  test  results,  the  slightly  above  aver- 
age (+)  and  well  above  average  (++)  performance  days  were 
combined  for  analysis  under  the  heading  of  above  average 
performance  days. 

a.   Test  for  Uniformity 

X2  Goodness -of -Fit  Test 

H  :   The  distribution  of  above  average  perform- 
ance  days  is  uniformly  distributed  over  the  intellectual 
biorhythm. 

H1 :   The  distribution  of  above  average  perform- 
ance days  is  not  uniformly  distributed  over  the  intellectual 
biorhythm. 

a  :   Level  of  Significance  =  .01. 


Interval  (days) 

Fr 

equency 

P. 

l 

nP. 

l 

1  <  t  <   5 

13 

5/33 

6.96 

5  <  t  <  10 

10 

5/33 

6.96 

10  <  t  <  15 

13 

5/33 

6.96 

15  <  t  <  20 

4 

5/33 

6.96 

20  <  t  <  25 

4 

5/33 

6.96 

25  <  t  <  33 

2 

8/33 

11.20 

X^(.99) 


15.  1 


V  =  21 


Table  X.   Intellectual  Data  (Above  Average  Performance  Days) 
X2  Goodness  -  of -Fit  Test  for  Uniformity. 


Decision:   Reject  H   in  favor  of  H, 
■"      o  1 
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Kolmogorov-Smirnov  Test 

Utilizing  the  same  H  ,  FL  ,  and  a  as  discussed 
&  o    1 


previously : 

Interval 

Frequency 

F0(x) 

S46^ 
|F  -  S| 


Kt<5  5<t<10  10<t<15  15<t<20  20<t<25  25<t<33 

13  10  13        4       4  2 

5/33  10/33  15/33  20/33  25/33  33/33 

13/46  23/46  36/46  40/46  44/46  46/46 

.1311  .1970  .3281  .2635  .3504  00.00 


Critical  D  =  .2403   Calculated  D  =  .3504 

Table  XI.   Intellectual  Data  (Above  Average  Performance  Days): 
Kolmogorov-Smirnov  Test  for  Uniformity. 

Decision:   Reject  H   in  favor  of  H1 . 

Both  of  the  above  tests  tested  for  uniformity. 

Lack  of  uniformity  suggests  the  presence  of  a  cyclic  phenomena, 

b.   Test  for  Distributional  Form 

X2  Goodness-of -Fit  Test 

H  :   The  distributional  form  of  above  average  per- 
o 

formance  days  is  N(8.5,5). 

hL  :   H   is  false . 
1    o 


a 


Level  of  Significance  =  .01 


Interval  (days) 

-  oo  <  t  <   3 

3  <  t  <  12 

12   <   t  <     CO 


Frequency 
7 

24 

15 


l 
1357 

6226 

2417 


nP. 
i 

6.24 


28.64 
11.12 


Table  XII.   Intellectual  Data  (Overall  Above  Average  Perform- 
ance Days):   x2  Goodness -of -Fit  Test  for  Distri- 
butional Form. 


X|C99)  =  9.21 


V  -  2.19 
Decision:   Do  not  reject  H 


4 .   Distribution  of  Below  Average  Performance  Days 

For  stronger  test  results  the  slightly  below  aver- 
age (-)  and  well  below  average  (--)  performance  days  were 
combined  for  analysis  under  the  heading  of  below  average 
performance  days. 

a.   Test  for  Uniformity 

X2  Goodness-of -Fit  Test 

H  :   The  distribution  of  below  average  perform- 
ance days  is  uniformly  distributed  over  the  intellectual 
biorhythm. 

H1 :   The  distribution  of  below  average  perform- 
ance days  is  not  uniformly  distributed  over  the  intellectual 
biorhythm. 

Level  of  Significance  =  .10. 


a 


Interval  (days)    Frequency   P. 


nP. 


1  <  t  <   8 

4 

8/33 

10.42 

8  <  t  <  16 

9 

8/33 

10.42 

16  <  t  <  27 

21 

11/33 

14.34 

27  <  t  <  33 

9 

6/33 

7.82 

X2(.90)  =  6.25 


V  =  7.42 


Table  XIII.   Intellectual  Data  (Below  Average  Performance  Days) 
X2  Goodness -of -Fit  Test  for  Uniformity. 

Decision:   Reject  PI   in  favor  of  H.  . 
J      o  1 
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Kolmogorov-Smirnov  Test 

Utilizing  the  same  H  ,  H1 ,  and  a  as  previously 


discussed 


Interval 

Frequency 

FG<X) 

S43OO 


o   43 


1  <  t  <  I 
4 
8/33 

4/43 

.1494 


8  <  t  <  16   16  <  t  <  27   27  <  t  <  33 
9  21  9 

16/33        27/33        33/33 


13/43 


.1825 


34/43 


.0275 


43/43 


00.00 


Critical  D  =  .1753   Calculated  D  =  .1825. 

Table  XIV..   Intellectual  Data  (Below  Average  Performance  Days) 
Kolmogorov-Smirnov  Test  for  Uniformity. 

Decision:   Reject  H  . 

J  o 

b.   Test  for  Distributional  Form 

X2  Goodness-of -Fit  Test 

H  :   The  distributional  form  of  below  average 
o  5 


performance  days  is  3(4,2) 


H-.  :H   is  false 
1   o 


a  :   Level  of  Significance 


Interval  (days)     Frequency 


01. 


P. 

1 


nP. 

1 


1  <  t  <  20 

13 

.350 

15.050 

20  <  t  <  30 

24 

.575 

24.725 

30  <  t  <  33 

6 

.075 

3.225 

X2(.99)  =  9.21      .V  =  2.69. 

Table  XV.   Intellectual  Data  (Overall  Below  Average  Perform- 
ance Days):   x?'  Goodness -of -Fit  Test  for  Distri- 
butional Form. 


Decision:   Do  not  reject  11 


B.   ANALYSIS  OF  ACCIDENT  DATA 

Sixty-six  accident  claims  from  the  Canadian  Forest  Prod- 
ucts Ltd.  (Howe  Sound  Pulp  Division)  were  analyzed  primarily 
with  respect  to  their  physical  biorhythms.   Those  accidents 
analyzed  were  single  victim  accidents  in  which  the  accident 
was  not  believed  to  be  caused  or  influenced  by  a  second 
party.   Table  XVI  summarizes  the  data  base. 

The  data  was  next  arranged  in  the  Physical  Biorhythm 
Data  Base  tableau  which  is  shown  as  Table  XVII. 

1.   Tests  for  the  Presence  of  a  Physical  Biorhythm 
a.   First  Analysis 

Since  accidents  are  considered  in  the  below 
average  performance  category,  the  x2  Contingency  Test  and 
X2  Coefficient  of  Contingency  are  inappropriate.   Conse- 
quently, the  analysis  begins  with  a  test  for  uniformity. 
X2  Goodness-of -Fit  Test 

H  :   The  occurrence  of  accidents  is  uniformly 
o  ' 

distributed  over  the  physical  biorhythm. 

H1 :   The  occurrence  of  accidents  is  not  uniformly 
distributed  over  the  physical  biorhythm. 

a  :   Level  of  Significance  =  .1. 


Interval  (days)     Frequency 


xJC90)  =  7.71 


V  =  8.63 


Table  XVIII 


Accident  Data  (1st  Analysis) 
Fit  Test  for  Uniformity. 


nP. 

i 


1  <  t  <  5 

10 

5/23 

14.35 

5  <  t  <  10 

08 

5/23 

14.35 

10  <  t  <  15 

20 

5/23 

14.35 

IS  <  t  <  20 

15 

5/23 

14.35 

20  <  t  <  23 

13 

3/23 

8.60 

X2  Goodness-of' 
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Subject 

Day 

of  Current    Day  of  Current 

Day  of  Current 

# 

Phys 

;  Biorhythm    Emot  Biorhythm 

Int   Biorhythm 

1 

21                22 

24 

2 

19                19 

26 

3 

13                22 

25 

4 

11                02 

18 

5 

03                02 

09 

6 

10                22 

14 

7 

21                21 

19 

8 

00                20 

17 

9 

19                11 

19 

10 

08                25 

19 

11 

19                11 

04 

12 

22                13 

18 

13 

20                05 

23 

14 

13                02 

21 

15 

07                11 

20 

16 

02                02 

10 

17 

04                24 

24 

18 

14                23 

21 

19 

16                24 

16 

20 

22                01 

07 

21 

16                15 

04 

22 

10                17 

07 

23 

11                18 

09 

24 

05                13 

21 

25 

05                09 

05 

26 

00                20 

17 

27 

14                02 

05 

28 

16                27 

25 

29 

21                27 

19 

30 

18                10 

24 

31 

19                09 

08 

32 

10                20 

05 

33 

20                19 

30 

34 

03                04 

14 

35 

22                10 

07 

36 

12                2  5 

06 

37 

00                01 

01 

38 

13                04 

21 

39 

21                20 

31 

40 

12                00 

27 

41 

17                02 

31 

42 

14                02 

11 

43 

21                10 

26 

44 

07                09 

08 

45 

11                19 

26 

Table  XVI.   Accident  Data  Base. 
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Table  XVI.   continued 


Subject  Day  of  Current  Day  of  Current  Day  of  Current 

#  Phys  Biorhythm  Emot  Biorhythm  Int   Biorhythm 

46  20  11  01 

47  14  17  24 

48  02  03  24 

49  21  09  17 

50  19  11  32 

51  01  15  16 

52  15  10  28 

53  17  09  04 

54  15  23  04 

55  00  18  24 


20 

11 

14 

17 

02 

03 

21 

09 

19 

11 

01 

15 

15 

10 

17 

09 

15 

23 

00 

18 

11 

15 

13 

01 

15 

11 

17 

26 

09 

14 

07 

15 

14 

25 

01 

16 

01 

09 

13 

18 

14 

04 

56  11  15  16 

57  13  01  07 

58  15  11  17 

59  17  26  06 

60  09  14  15 

61  07  15  12 

62  14  25  15 

63  01  16  07 

64  01  09  04 

65  13  18  51 

66  14  04  20 
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Cycle  Day 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

00 

Accident  Frequency 


3 
2 
2 

1 
2 

3 
1 
1 
3 
4 


2 
5 
6 
3 
3 
3 
1 
5 
3 
6 
3 
4 


Frequency  refers  to  the  total  number  of  accidents  which 
occurred  on  the  respective  biorhythm  day.   For  example,  of 
the  66  subjects,  three  experienced  accidents  on  the  first 
day  of  their  physical  biorhythm,  two  experienced  their  ac- 
cidents on  the  second  day,  etc.   The  data  can  now  be  easily 
categorized  for  use  in  the  various  statistical  models. 


Table  XVII.   Physical  Biorhythm  Data  Base 
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Decision:   Reject  H   in  favor  of  H, . 
J      o  1 

Kolmogorov-Smirnov  Test 

Utilizing  the  same  H  ,  KL  ,  and  a  as  previously 


discussed 


Interval  (days) 

Frequency 

F0(x) 

S66^ 


o   66 


l<t<5  5<t<10  10<t<15  15<t<20  20<t<23 

10  8  20  15  13 

5/23  10/23  15/23  -20/23  23/23 

10/66  18/66  38/66  53/66  66/66 

.0655  .1620  .0764  .0665  00.00 


Critical  D  =  .1502 


Calculated  D  =  .1620. 


Table  XIX.   Accident  Data  (1st  Analysis) :   Kolmogorov-Smirnov 
Test  for  Uniformity. 

Decision:   Reject  H   in  favor  of  H1 . 
J  o  1 

Rejection  of  the  uniformity  hypothesis  by  both 
statistical  tests  suggests  the  presence  of  a  cyclic  phenomena, 
b.   Second  Analysis 

The  Theory  of  Biorhythms  has  a  built  in  variance 
of  ±12  hours  (24  hours)  due  to  the  uncertainty  of  the  exact 
hour  of  birth.   Consequently,  it  is  possible  that  the  original 
calculations  may  be  as  much  as  24  hours  away  from  where  they 
appear  in  the  original  data  tableau.   This  analysis  investi- 
gates the  statistical  significance  in  shifting  the  data  24 
hours  (optimally) . 

X2  Goodness-of -Fit  Test 

H  :   The  occurrence  of  accidents  is  uniformly 
o 

distributed  over  the  physical  biorhythm. 
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Cycle  Day 

Accident  Frequency 

1 

1 

2 

1 

3 

2 

4 

1 

5 

1 

6 

1 

7 

2 

8 

2 

9 

1 

10 

1 

11 

2 

12 

6 

13 

5 

14 

6 

15 

3 

16 

3 

17 

3 

18 

1 

19 

5 

20 

3 

21 

6 

22 

3 

23 

7 

Table  XX.   Optimal  Data  Tableau 


H 


. :   The  occurrence  of  accidents  is  not  uniformly 


distributed  over  the  physical  biorhythm. 

a  :   Level  of  Significance  =  .005. 

Interval  (days)     Frequency    P. 


X|C905) 


14.9 


19.71. 


nP. 

l 


1  <  t  <  5 

6 

5/23 

14.35 

5  <  t  <  10 

7 

5/23 

14.35 

10  <  t  <  15 

22 

5/23 

14.35 

35  <  t  <  20 

15 

5/23 

14.  35 

20  <  t  <  23 

16 

3/23 

8.60 

Table  XXI.   Accident  Data  (2nd  Analysis) 
Fit  Test  for  Uniformity. 


X2  Goodness-of- 
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Decision:   Reject  H   in  favor  of  H..  . 
J      o  1 

Kolmogorov-Smirnov  Test 

Utilizing  the  same  H  ,  H...  ,  and  a  as  previously 


d 

iscussed 

Interval 

Frequency 

F0Cx) 

S66^ 

F   -  S££ 

66 


l<t<5   5<t<10  10<t<15  15<t<20  20<t<23 

6       7  22  15  16 

5/23   10/23  15/23  20/23  23/23 

6/66   13/66  35/66  50/66  66/66 

.1265  .2378  .1219  .1120  00.00 


Critical  D  =  .2006    Calculated  D 


2378 


Table  XXII.   Accident  Data  (2nd  Analysis);-  Kolmogorov-Smirnov 
Test  for  Uniformity. 

Decision:   Reject  H  . 
J      o 

The  above  analysis  shows  that  shifting  the  oc- 
currence within  the  cycle  by  24  hours  yields  more  significant 
results.   However,  knowledge  of  the  subject's  exact  hour  of 
birth  eliminates  the  need  for  shifting  the  data  24  hours  and 
yields  more  reliable  results. 
c.   Third  Analysis 

This  analysis  dichotomizes  the  data  into  loss- 
time  accidents  and  non-loss -time  accidents.   A  loss-time 
accident  is  an  accident  which  results  in  the  employee  re- 
maining out  of  work  for  at  least  one  day.   The  following 
tableau  contains  data  from  loss -time-accident  subjects. 
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Subject 


01 
11 
18 
20 
21 
23 
39 
41 
53 
58 
66 


Day  of  Current 
Phys  Biorhythm 

*21 

19 

14 
*22 

16 
*11 
*21 

17 

17 

15 

14 


Day  of  Current 
Emot  Biorhythm 

22 

11 

23 
*01 
*15 

18 

20 
*02 

09 

11 

04 


Day  of  Current 
Int   Biorhythm 

24 

04 

21 

07 

04 

09 
*31 
*31 

04 
*17 

20 


*It  is  possible,  assuming  ±12  hours  variance, 
that  these  days  could  have  been  critical  days. 

Table  XXIII.   Loss-Time  Accidents. 

X2  Contingency  Test 

H  :   Loss-time  accidents  and  non- loss -time  acci- 
o 

dents  independent  of  where  they  occur  within  the  physical 
biorhythm. 

H1 :   Loss-time  accidents  and  non- loss -time  ac- 
cidents are  dependent  of  where  they  occur  within  the  physical 
biorhythm. 

a  :   Level  of  Significance  =  .05. 


Loss -Time 
Accidents 


Non-Loss -Time 
Accidents 


1st  h   Cycle 
2nd  h.   Cycle 


1 

21 

11 

33 

X2(.95)  =  3.84 


V  -  4.18 


Table  XXIV.   x2  Contingency  Test. 

Decision:   Reject  H   in  favor  of  H. 
J  o  1 


6  3 


This  result  illustrates  more  severe  accidents 
are  dependent  upon  the  sequential  day  of  the  subjects' 
biorhythm. 

X2  Coefficient  of  Contingency 


C-       ■     Y 


\V   +  N 


r  -  1    4.18 
u     \J4. 18  +  66" 

C  =    .244. 

This  represents  the  degree  of  association  or  the 
strength  of  the  relationship  previously  measured  by  the  x2 
Contingency  Test.   The  Coefficient  of  Contingency  is  just 
another  way  of  examining  the  strength  of  the  rejection  in 
the  x2  Contingency  Test.   Had  the  x2  Contingency's  null 
hypothesis  (that  of  independence)  has  been  rejected,  there 
would  have  been  no  reason  to  calculate  a  Contingency  Coeffi- 
cient since  the  x2  deviation  statistic  would  have  been  so 
small  that  the  resulting  Coefficient  would  have  been  insig- 
nificant.  This  can  be  illustrated  by  examining  the  formula. 
In  order  not  to  reject  the  null  hypothesis,  the  deviation 
statistic  would  have  to  be  less  than  3.84.   Assume  that  this 
statistic  was  1.5.   With  the  same  sample  size  of  66,  the 
Coefficient  of  Contingency  is  equal  to  .149  which  is  insig- 
nificant at  any  alpha  level.   A  close  examination  of  the 
parameters  V  and  N  illustrates  that  one  must  be  careful  when 
interpreting  this  coefficient.   A  plot  of  the  Contingency 
Coefficient  against  sample  size,  holding  the  deviation  constant 
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will  show  that  as  the  sample  size  becomes  large,  the  coeffi- 
cient becomes  asymptotic  to  zero.   Consequently,  it  is  more 
sensitive  to  smaller  sample  sizes,  as  is  the  case  with  this 
biorhythm  data. 

2.   Criticality 

Lack  of  data  eliminates  an  indepth  analysis  of 
criticality,  however,  it  is  interesting  to  observe  the 
occurrence  of  criticalities  among  loss- time  accident  victims. 
Let 

X  =  criticalities 
P(PCD)  =  .1739 
P(ECD)  =  .1428 
P(ICD)  =  .1212. 


Biorhythm 

E(X) 

V(X) 

#  Critical 

Observations 

Physical 

1.91 

1.58 

4 

Emotional 

1.57 

1.35 

3 

Intellectual 

1.33 

1.16 

3 

Table  XXV.   Accident  Data  Criticality  Tableau. 

The  following  probabilities  were  calculated  using 
the  binomial  model. 
Let 

X  =  #  criticalities  observed  within  a  physical 
biorhythm  out  of  11  possible  observations. 
P(PCD)  -  .1739. 
Table  XXVI  illustrates  that  the  expected  (mean) 
number  of  observed  criticalities  (1.91)  occurs  with  highest 
probability  while  four  or  more  criticalities  can  be  considered 
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a  rare  event.  Figure  5  illustrates  that  90%  of  the  proba- 
bility of  occurrence  of  critical  observations  is  contained 
within  1.5  standard  deviations  of  the  mean. 

n      P(X=n) 

0  .123300000 

1  .283109200 

2  .297651200 

3  .187889600 

4  .077992200 

5  .023344900 

6  .004905300 

7  .000722300 

8  .000074400 

9  .000003700 

10  .000000158 

11  .000000003 

Table  XXVI.   Criticality  Prob  Tableau. 
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9   10    11 


Vi- 1.50  -'y 


y+1.50 


CRITICALITIES 


Figure  5.   Observed  Accident  Data  Density  Function 
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VI.   DISCUSSION 

Results  based  on  statistical  models  are  heavily  dependent 
on  the  researcher's  interpretation  of  the  theory.   One  who 
considers  criticality  as  the  backbone  of  the  theory  will 
stress  different  aspects  of  biorhythmic  analysis  than  will 
a  researcher  who  believes  criticality  plays  a  minor  role  in 
biorhythmic  prediction.   Consequently,  the  choice  and  de- 
velopment of  the  combined  statistical  model  illustrates  that 
every  aspect  of  a  biorhythm  should  be  analyzed  prior  to  draw- 
ing conclusions.   Use  of  the  models  in  a  thorough  and  rigid 
analysis  would  allow  one  to  draw  objective  conclusions  based 
on  fact  instead  of  formulating  broad  generalizations  based 
on  subjective  interpretations  of  the  theory. 

The  analysis  of  intellectual  and  accident  data  are  not 
meant  to  validate  the  theory.   The  analyses  both  illustrate 
the  use  of  the  statistical  models  as  well  as  support  this 
paper's  interpretation  of  biorhythms.   The  analysis  of  in- 
tellectual data  rigidly  adhered  to  the  combined  statistical 
model  for  testing  the  presence  of  a  biorhythm.   At  the  .005 
level  of  significance  the  Chi-Square  Contingency  Test  es- 
tablished that  a  relationship  existed  between  the  different 
levels  of  performance  and  the  days  on  which  they  occurred. 
Supporting  this  result  was  an  impressive  .48  coefficient  of 
contingency.   As  was  expected  the  application  of  both  the 
Chi-Square  Goodness -of -Fit  Test  and  the  Kolmogorov- Smi rnov 
Test  (KS)  to  average  performance  data  showed  that  days  on 
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which  the  subject  was  not  intellectually  challenged  (average 
performance  days)  were  uniformly  distributed  over  the  entire 
cycle.   This  result  is  coincidental  with  biorhythm  theory. 
Biorhythms  measure  variation  in  performance  potential.   On 
average  days  variation  is  zero  and  thus  unpredictable  and 
completely  random  as  in  the  case  of  the  uniform  distribution, 
However,  repeated  application  of  the  Chi-Square  and  KS  Tests 
to  the  above  and  below  average  performance  days  illustrated 
non-uniformity  which,  according  to  the  model  implies  a 
cyclic  phenomena.   Identification  of  the  distribution  of 
both  the  above  and  below  average  performance  days  was  pos- 
sible.  Identifying  the  distribution  is  instrumental  in  that 
it  allows  one  to  determine  expected  frequencies  and  their 
variations.   Although  application  of  both  statistical  tests 
was  not  necessary,  the  use  of  the  additional  test  serves  as 
a  check  and  lends  supporting  evidence  to  the  original  analy- 
sis.  For  small  samples  the  KS  Test  is  more  powerful  since 
the  Chi-Square  Test  is  forced  to  combine  categories  when 
expected  value  criteria  are  not  met.   It  is  realized  that 
the  limited  sample  size  presents  difficulties  when  attempt- 
ing to  dichotomize  in  attempts  to  gain  as  much  information 
as  possible  from  the  data.   However  the  significance  of  the 
results  is  impressive  regardless  of  the  sample  size.   The 
presence  of  only  one  critical  event  prevented  a  criticality 
analysis  . 

There  are  many  problems  encountered  when  attempting  to 
analyze  data  such  as  academic  grades.   It  lias  been   the 
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author's  experience  that  exams  on  the  graduate  level  are 
designed  to  challenge  the  student.   Most  conscienious  stu- 
dents over-prepare  for  exams  thus  reducing  the  exams'  chal- 
lenge.  This  makes  it  difficult  to  collect  and  analyze  clean 
data.   If  the  intellectual  biorhythm  actually  measures 
memory  potential  then  a  knowledge  of  potential  on  exam  days 
may  be  a  determinant  in  how  much  preparation  is.  required 
for  the  exam. 

In  general,  analyses  of  accident  data  is  not  considered 
appropriate  in  analyzing  biorhythms  since  accident  data 
isolates  the  analysis  to  the  low  performance  potential  area 
of  the  biorhythm.   This  single  class  of  data  prevents  the 
use  of  the  entire  statistical  model.   For  example  the  Chi- 
Square  Contingency  Test  and  the  Coefficient  of  Contingency 
are  not  applicable  as  suggested.   However  application  of 
both  the  Chi-Square  Goodness -of -Fit  Test  and  the  KS  Test 
established  the  existence  of  a  non-uniform  distribution  of 
below  average  performance  days.   One  may  wonder  what  is  the 
effect  of  considering  a  natural  built  in  variance  due  to 
the  uncertainty  of  the  exact  hour  of  birth  and  possibly  the 
exact  time  of  the  observed  event.   Between  critical  days 
(above  and  below  average  performance)  the  variance  has  a 
negligible  effect  on  whether  the  biorhythm  is  present.   How- 
ever it  has  a  profound  affect  on  the  criticality  in  that  the 
assumption  of  variance  increases  the  chance  probability  of 
an  event  occurring  on  a  critical  day.   Further  dichotomizing 
the  time-loss  and  non-time- loss  accidents  allowed  the  assessment 
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of  the  severity  of  the  accidents.   The  refinement  in  classi- 
fication resulted  in  applying  the  Chi-Square  Contingency 
Test  and  the  Coefficient  of  Contingency.   This  showed  that 
time-loss  accidents  (more  severe)  were  dependent  upon  the 
biorhythmic  day  on  which  they  occurred.   It  is  also  sug- 
gested that  the  occurrence  of  critical  events  is  signifi- 
cant.  For  a  sample  size  of  11  the  expected  number  of  events 
occurring  on  a  critical  day  is  1.91  with  a  variance  of  1.58 
and  a  standard  deviation  of  1.25.   The  fact  that  four  criti- 
cal observations,  lie    outside  1.5  standard  deviations  of 
the  mean  suggests  the  presence  of  the  criticality  phenomena. 
Although  it  was  possible  to  calculate  the  exact  probability 
of  this  event,  it  is  felt  that  further  attempts  at  accuracy 
will  not  be  meaningful. 

P(3  phy  critical  days;  1  double  critical  day; 
7  non-critical  phy  days)  = 

11  3  7 

.  _   (.1739)-5(.0124)  (.8137)7  =  .02. 

As  illustrated,  this  can  be  considered  a  rare  event.   On  the 
average  it  is  expected  that  1.91  events  would  occur  on  criti 
cal  days.   It  is  also  expected  that,  under  a  null  hypothesis 
of  pure  chance,  the  expected  number  is  not  rigid  and  has  a 
variance  of  1.58  and  a  standard  deviation  of  1.25.   This 
means  that  if  the  occurrence  of  critical  events  was  pure 
chance  one  could  expect  as  few  as  zero  observed  critical 
events  or  as  many  as  3.82.   The  fact  that  four  critical 
events  were  observed  suggests  that  criticality  does  not 
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occur  in  a  pure  chance  fashion  since  more  critical  events 
occurred  than  were  expected.   Consequently  it  is  felt  that 
a  criticality  phenomena  is  present. 

Although  both  analyses  do  not  comprise  a  validation, 
each  analysis  completely  supports  the  theory  in  a  very 
strong  sense. 
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VII.   CONCLUSION 

Although  biorhythmic  phenomena  is  considered  esoteric 
and  attributed  to  the  occult  sciences,  the  possibility  that 
such  isochronic  phenomena  exists  should  not  be  dismissed. 
As  suggested  there  are  scientific  methodological  ways  of 
investigating  the  theory.   The  method  presented  here  is  not 
considered  the  only  method.   Preliminary  investigations 
suggest  the  validity  of  the  theory,  however,  only  a  scien- 
tifically based  validation  will  be  accepted  by  a  scientific 
community . 
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